Rolling friction of a hard cylinder on a viscous plane 
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The resistance against rolling of a rigid cylinder on a flat viscous surface is investigated. We found 
that the rolling-friction coefficient reveals strongly non-linear dependence on the cylinder's velocity. 
For low velocity the rolling-friction coefficient rises with velocity due to increasing deformation rate 
of the surface. For larger velocity, however, it decreases with velocity according to decreasing contact 
area and deformation of the surface. 

PACS numbers: 46.30.Pa, 62.40. +i, 8f.40.Pq 
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I. INTRODUCTION 

The effect of rolling friction has been investigated by 
many scientists according to its great importance in en- 
gineering (e.g. Q) and physical science (e.g. jl). Sci- 
entific publications on rolling friction range back to, at 
least, 1785 when Vince described systematic experiments 
to determine the nature of friction laws ||. 

It is known that surface effects such as adhesion 
(e.g. (J]), electrostatic interaction (e.g. ||), and other 
surface properties (e.g. 0]) may have strong influence 
on rolling friction. For viscoelastic materials, however, 
it was argued that rolling friction is due very little to 
surface interactions, i.e. the major part is due to de- 
formation losses within the bulk of the material 
Under this assumption Greenwood et al. calculated 
the rolling friction coefficient for a hard sphere rolling on 
a soft plane. The deformation in the bulk was assumed 
to be completely plastic. Then an empirical coefficient 
was introduced to account for the incomplete recover of 
the material. Recently a similar problem has been ad- 
dressed in H where the rolling friction coefficient for a 
soft sphere on a hard plane has been derived as a first- 
principle continuum-mechanics expression. This coeffi- 
cient has been found within a quasi-static approach |l(|] 
as a function of the viscous and elastic constants of the 
sphere material without introducing phenomenological 
parameters. 

In the case of a soft sphere rolling on a hard plane ||] 
the contact surface between the bodies is flat. This allows 
for the application of Hertz's contact theory. In the oppo- 
site case of a hard sphere or cylinder on a viscous plane 
which we address here, this assumption is not justified 
since the plane deforms in such a way that its shape fol- 
lows the shape of the rolling body in the area of contact. 
This complicated shape of the contact surface excludes 
the direct application of Hertz's contact theory and may 
violate Hertz's contact law which relates the force acting 
between the interacting bodies to their deformation. 

The velocity dependence of the rolling friction coeffi- 
cient originates from the fact that the deformation of the 
surface varies with the velocity of the rolling body. For 



small velocities the viscous stress, proportional to the de- 
formation rate, is small. In this case the deformation of 
the plane (measured by a depth h of which the body pene- 
trates the surface, see Fig. |l|) is determined mainly by the 
elasticity of the plane and by the weight of the body. On 
the other hand at very large velocities the viscous stress 
becomes comparable to the elastic stress. As a result the 
plane supports the rolling body at significantly smaller 
deformations. This leads to decreasing penetration depth 
h, hence, less energy may be required to deform the sur- 
face. In this case one observes decreasing resistance to 
rolling with increasing velocity. Mainly because of the 
complicatedly shaped contact area it is not possible to 
treat these effects within the first-principle continuum- 
mechanics description, hence a simpler model will be con- 
sidered. As shown in the following this model reflects 
the most important properties of the problem addressed. 
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FIG. 1. Sketch of the system: A rigid cylinder rolls 
on a plane built up of independently moving damped 
springs. According to the motion of the cylinder the sur- 
face is deformed and mechanical energy is lost due to 
its damped motion. This dissipation of energy can be 
understood as rolling friction. 
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II. THE MODEL 

We investigate the resistance against rolling of a hard 
cylinder on a soft plane. Consider a cylinder of radius R, 
length L 3> R and mass per unit length M which rolls 
along the x-axis with velocity v. We assume that the 
surface may be modelled by a system of noninteracting 
springs. Their elastic, damping and inertial properties 
are described by the coefficients k, 7 and m. These are 
defined as follows: kdx and jdx give the elastic and 
viscous force per unit length along the cylinder axis, 
m dx gives correspondingly the linear mass density of 
the springs (i.e. their mass per unit length along this 
axis). The viscous plane becomes deformed in the range 
x- < x < x + due to the mass of the cylinder. For defor- 
mation rates small as compared with the vertical speed 
of sound, Hooke's law is valid and one has the equation 
of motion which describes the behaviour of the surface 



my{x) + -yy(x) + ky(x) = f(x,t) , 



(1) 



where / is the force density (/ dx gives the force per unit 
length along the cylinder axis) which acts on the plane 
in the region of contact. Outside of the contact area 
the force density is zero. In our model we neglect lateral 
couplings of the springs (for justification of these assump- 
tions see the Appendix). We also assume that no surface 
waves are excited on the plane, i.e. that the condition of 
the overdamped motion of the surface 



m 1 7 
~k < 4 P 



(2) 



holds true [0|. Since we are interested in the steady 
state, i.e. time independent behavior, we do not need 
to explicitly consider the tangential interactions between 
the cylinder and the surface. An arbitrarily small tan- 
gential force assures that the cylinder does not slide. If 
one, however, is interested in accelerated motion of the 
cylinder, tangential forces need to be considered. 

By geometrical considerations we find for the shape 
y(x) of the surface of the deformed plane in the contact 
area: 



y(x) = R— h — \J R? — (x — x c ) 2 , X- < x < x + , (3) 

where x c is the x-coordinate of the center and h — 
—y-min = —y{x c ) is the penetration depth of the cylin- 
der. For \x - 



<C R we approximate (||) by 



(4) 



The center of the cylinder moves with constant velocity 
v, i.e. x c = vt. Hence, the time derivatives of y{x) read 



y(x) = -x, 



R 



R 



y(x) = -j^x c = — = const. 



(5) 
(6) 



The compression force exerted by the plane to the 
cylinder is 



x + 



F, 



N 



f(x)dx. 



(7) 



(For simplicity of notation here and in what follows we 
notate the forces, energy and torque per unit length of 
the cylinder, e.g. the total force exerted by the plane to 
the cylinder is LFjy.) 

The springs at 2+ which at time t get in contact with 
the cylinder need a separate discussion |L2| : At time t — S 
(5 — > 0) their velocity is zero while infinitesimal time later 
y(t + 5) is finite according to (S). This singularity in the 
velocity distribution may be attributed to a force F' N , 
acting at point x+. This gives a finite contribution to 
the total force which can be determined by the following 
consideration: 

In the time interval dt the cylinder moves by vdt. So 
it accelerates springs of total mass mvdt. The total mo- 
mentum received by these springs is dp — y(x + ) mvdt, 
hence, 

F N = --^ = -y(x + )mv = — mv . (8) 

The total force Fn + F^ supports thus the weight of the 
cylinder: 

F N + F' N = Mg . (9) 

Substituting y and its time derivatives (Eqs. (|],[|]6|)) 
in (Q) we get an expression for the force density in the 
contact area 

f( x , t) = A {x - Xc f -ll^-^ + ^-hk ( io) 

which has to satisfy the contact condition 

f(x,t)<0 for x_ < x < x + (11) 

and which determines the boundaries of the contact area 
x± . In comoving coordinates £ = x — x c Eq. (jig) reads 



J ^ J 2iT R^ R 



(12) 



The boundary of the contact area at the front side of the 
cylinder in the direction of motion is 



£+ = V2Kh 



(13) 



according to geometry. The boundary at the back side is 
determined by /(£-) — 0, i.e. 



k 2 iv mv 2 

0= 2Rt---Rt- + ^r~ hk 



(14) 
(15) 
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For the self-consistency one needs the condition £_ < 
V2Rh, i.e. h(v) > m 2 v 2 /2^ 2 R to be fulfilled, which re- 
stricts the velocity v from above 

Because of the dissipative properties of the plane the 
motion of the cylinder corresponds to a loss of mechanical 
energy per time. Another contribution originates from 
the instantaneous acceleration of plane material which 
gets in contact with the plane during the time dt p2[ . 
Therefore, we find for the energy which is transferred 
per time from the cylinder to the plane 

E = - ( dtfiom - m^l v 



Substituting Eqs. (|l|,|o|) into Eq. (|T|) and using (|J|) one 
can solve it perturbatively to find the front boundary 



where 



3R \ 1/3 Xv 2 



4£ 



(0) 



(21) 



1/3 



(22) 



denotes the front boundary for the static case. Then from 
(p0|) the rear boundary follows 



|/^/(0-^ 3 | 

i- 
-vF R . 



(16) 



Equation (|16| ) defines the force Fr which acts against 
rolling of the cylinder 



F R = ~Jdtim)+rnv 2 ^. 



(17) 



£ --£V>+v?-- — + ->> 
4 - _ § + +V k 4^(0) + ■ 

Correspondingly, the penetration depth reads: 



where 



h Q ^h^ = [^) 2 /2R 



(23) 



(24) 



(25) 



To evaluate Fr we need the force /(£) given by Eq.jl^) 
with the penetration depth h(y) which is to be found. 
The penetration h(v) results from equilibrating the dy- 
namical resistance of the surface with the weight of the 
cylinder (^|) 



Mg = - J 



2R 



e- 



71) 



- hk 



R 



-V2hR 



(18) 



is an 



with the boundaries ( p"3| ) and (|T^). Equation 
implicit equation for h. 

We consider first the limit of small velocities. In this 
limit Eq. (|18| ) may be solved as a perturbation expansion, 
h = +h> 1 ^v + h^v 2 +••: It is more convenient, how- 
ever, to solve (|l^) with respect to an expansion of £ + : 



(0) 



v^+v 2 ^ 



(19) 



Using Eqs.(|l3|), (|T|) and © one can further write the 
small- v expansion for : 



(20) 



where 



is the static penetration depth. From Eqs. (J2j) and (g) 
it follows that the penetration depth h decreases with 
increasing velocity. 

Using the obtained expansions for £ + and £_ it is 
straightforward to calculate the rolling friction force. 
Substituting (|Io|), (|l]) and (||) into Eq. (0) one finally 
arrives at an expression for the rolling friction torque, 
M = RF R : 



M = UrouMg 



Mroll 
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3A 
4^f 



(26) 
(27) 



As it follows from Eq. (^7|) , in the limit of small veloc- 
ities the leading linear term depends only on the viscous 
and elastic constants and does not depend on the iner- 
tial properties of the material, characterized by m. This 
means that in this regime the inertial effects in the de- 
formation process of the plane are negligible. The second 
nonlinear term takes into account via m/k in A the iner- 
tial effects up to 0{v 2 ). It is also interesting to note, that 
while the linear term does not depend on the radius of 
the body, the nonlinear term depends (via ) on both, 
the radius of the cylinder and its mass per unit length 
M. 

In the general case Eq. ( |l8| ) has to be solved numeri- 
cally. The velocity dependence of the penetration depth 
h and of the rolling friction coefficient ^ TO \\ are shown in 
Fig. | and Fig. [| Calculations were performed for a steel 
cylinder of radius R = 0.1 m, and of mass per unit length 
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M = 250 kg m 1 rolling on the rubber surface with 
the following material parameters: m = 100 kg • m , 



= 10 7 kg • m 2 • s 2 , 7 = 5 ■ 10 5 kg • m 



For velocities v > v* the rolling friction coefficient de- 
creases with increasing velocity. For particular param- 



These eters used here one obtains \x 



values were obtained from density and elastic constants of 
rubber |Q ; the viscous constant was estimated from the 
restitution coefficient for colliding rubber spheres, simi- 
larly as in H, where the details of this estimate are given. 
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FIG. 2. The penetration depth ft as a function of the 
velocity v according to the numerical solution of Eq. ( |l8| . 
The dashed line shows the approximation (|24[). 
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v [m/sec] 

FIG. 3. The rolling friction coefficient /i ro n over the 
velocity v (numerical solution). The dashed line shows 
the approximation (p7[). 

As shown in the figures, at high velocities the depth 
h decreases and the rolling friction coefficient reaches a 
maximal value of 



A'roll "n u 



3 k 2 X 



(<>) 



at 



27^ 
3 k X 



(28) 



(29) 



roll 



1.125 ■ 10" 2 m and 
v* = 0.4500 ms -1 , which are in a reasonable agreement 
with the numerically obtained values (see Figs, g and|^). 

Note that with increasing velocity the rear boundary 
£_ shifts in positive direction and approaches the front 
boundary £ + . The contact area between the surface and 
the cylinder thus decreases and at some critical velocity 
it shrinks to a point (a line along the cylinder axis). For 
velocities larger than this critical one the weight of the 
cylinder is sustained by the force F' N , acting at a single 
point at the front boundary £+. 



III. SUMMARY AND DISCUSSION 

We investigated the rolling motion of a hard cylinder 
on a viscous plane. The elastic, viscous and inertial prop- 
erties of the plane were modelled by a system of uncou- 
pled springs which are characterized by linear elastic and 
viscous coefficients and mass density. For small velocity 
of the rolling cylinder we determined the velocity expan- 
sion of the rolling friction coefficient analytically up to 
second order. For larger velocities the rolling friction 
was determined numerically. 

In the range of low velocity our analysis shows increas- 
ing rolling friction coefficient with increasing velocity. At 
a certain velocity the coefficient reaches its maximal value 
and decreases when the velocity is further increased. 

In the low velocity regime where the rolling friction 
coefficient rises linearly with the velocity of the cylinder 
its value depends on the ratio of the viscous and elastic 
constants. For high velocities one notes nonlinear depen- 
dence on the cylinder's mass, radius and on the material 
constants of the surface. 

We also analyzed the dependence of the penetration 
depth (i.e. the depth at which the cylinder sinks into 
the bulk of the surface) on the cylinder velocity. We 
found that the penetration depth decreases with increas- 
ing velocity, i.e. that the rolling cylinder emerges up- 
wards when its velocity increases. 
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APPENDIX: LATERAL COUPLINGS 

In our model we assumed that the viscous surface is 
composed of springs which are not coupled in lateral (hor- 
izontal) direction. This simplification may limit the va- 
lidity of the model since we ignore lateral interactions 
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in the plane, which obviously exist in real systems and 
which cause the tangential elasticity of the surface. In 
the appendix we want to study in more detail the range 
of validity of our model. 

Lateral interactions of the springs may be taken into 
account if instead of Eq. ([!]) one uses 

d 2 

my(x) + -yy(x) + ky(x) - d-^y = f{x, t) , (Al) 

where the constant d describes the lateral coupling be- 
tween the springs. With the dimensionless variables 
y = y/R , x — x/R and i —tj (R/v) the left-hand side of 
Eq. (Al) reads 



ttiv 2 2 d 



RdZ 



2 y(x). (A2) 



Thus we conclude that if the radius of the cylinder sat- 
isfies the condition 



R 2 » d j 
k 



(A3) 



the term describing the lateral coupling is negligible and 
our model is valid. 

The value of \fdfk is a characteristic length of influ- 
ence of lateral couplings. If the above condition (A3) 
is satisfied, one can apply Eq. ([!]) for almost the entire 
contact region, except for the (small) region around £+. 
In these region the deformation of the surface may differ 
from the idealized shape which, as we assumed, results 
from pure geometrical considerations. Instead the range 
of deformation of the surface is slightly enhanced by a 
characteristic size ~ rg due to lateral couplings between 
the springs. If the cylinder is at rest or if it moves very 
slowly, similar discussion applies to the small region ~ r 
around the rear point . Consider first the motion of 
the surface in a region around £+. 

Since tq the region of the characteristic size ~ ro 

may be considered as a point, so that the total force 
acting on this region may be attributed to the single 
point £+. The dissipation in this region (which is finite) 
is, therefore, attributed to the point £ + too. Coarse- 
graining, therefore, results in the force (£ + ) acting at 
the point £+. 

We will estimate the characteristic size ro: If we take 
into account lateral couplings in the region ro around £ + 
there is no contact between the surface and the cylin- 
der (free surface condition). The derivative y' = dy/dx 
changes from zero (condition of the undisturbed plane) to 
v' ~ ?/(£+) ~ £+/R a t the point of contact, which follows 
from the geometry of the system (see Fig. [|). Thus, in 
this region one estimates y" = d 2 y/dx 2 ~ £ + /i?ro- Sim- 
ilarly, the characteristic value of y in this region reads 
y ~ r o • y' ~ r o ■ (,+ / R and, correspondingly, the charac- 
teristic values y = y(£+) = —v^ + /R and y = v 2 / R follow 
from Eqs. (§), (§). Then we write the condition of the 
free surface in this transient region, 



mv 2 /R - 7<+/i? + kr £ + /R ~ d ■ £+/Rr = , (A4) 



to estimate the size of the region: 




/ mv 2 




f mv 2 






2k J 


\2k£+ 


2k J 



(A5) 



which yields r = Wd/k for the static case. Hence, the 
condition for coarse-graining, ro <C £+, reads 



m 



(A6) 



We want to discuss the consequences of the assumption 
QA6| ) or of the assumption ^ ^> d/k (the later condition 
follows from the former o ne, u nless the velocity is too 
high, i.e. unless v ^> £+y/k/m). In this case the force 
Fjf(£+) which acts at the point £ + reads 



(A7) 



5+ -5 



where 5 is of the order (say somewhat larger) than ro, 
and we can wri te for the different terms in the left-hand 
side of Eq. Q: 



mydx 



mv 



i++s d y, 

—ax 

i+-s dt 



-+S 



, dx 



-5 dt 



e++<5 



dy = mv [y(£ + + S) - y(£ + ~ 5)] 



-mvy(£+) = mv 2 £+/R 



(A8) 



where we take into account that y(£,+ + 5) = (the sur- 
face is at rest) and that y(£+ — S) = y{£+) on the coarse- 
grained scale. Similarly, using the above estimate of y(x) 
in the transient region, we obtain the coarse-grained es- 
timates: 



r6++S 

7 / ydx = ~/v [y(£ + + 5) - y(£ + - 5)] 



-Jwr Q £ + /R 
(A9) 



and 

/•£++« 
k I ydx ' 
J €+-6 



k y'(£ + )xdx ~ 2kr^\/R. (A10) 

h+-6 



Finally, the last term reads 



-d- I y"dx 
--d-y'^ + ) = d-i + /R. 



d-[y'(Z + +S)-yi(Z + -6)} 

(All) 



As it follows from Eqs. ( A8|- All) the second and third 
terms, proportional to ro, vanish on the coarse-grained 
scale. The fourth term does contribute to Fj^(£+) on the 
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coarse-grained level, but it does not depend on the ve- 
locity v. It may be taken into account within the general 
scheme of calculation of the rolling friction given above. 
Namely, with this term included, one obtains, e.g. for 
the front boundary 



(0) _ i 3R Mg 

4+ 2k 



1 + 77 



(d/k) 



1/3 



^ 3R M g x 1/3 



2/v 



(A12) 



where the condition ^> d/k was used. Similarly, 
the impact of this term on the other expressions ob- 
tained previously is negligible, i.e. it is of the order 

(d/k)/ Oz+j ^ 1 under the coarse-grained condition. 

Thus, we conclude that the fourth term in F' N (£+), which 
accounts for the lateral interactions may be also ne- 
glected. This gives the result 



F' N ^ + )=mv 2 i+/R. 



(A13) 



Similar considerations may be performed for the dis- 
sipation in the transient region. Skipping the details of 
the analysis (very similar to that for F' N (£ + )), we give 
the final result: 



y 2 (M 



(A14) 



which describes the energy loss in the point £ + on the 
coarse-grained scale. 

Using the same reasoning one can consider the region 
~ ro around the rear point £_ to conclude that under the 
condition 3> d/k its contribution to the total force and 
to the dissipation is negligible: Indeed, the region ~ rg in 
the rear part of the contact area may influence the mo- 
tion of the cylinder only for very small velocities v when 
| — |£_| ~ tq. However, contribution to the force and 
dissipation from this rear region is proportional to ro and 
therefore may be neglected. 

Thus, we conclude that our simplified model of the vis- 
cous surface as a system of linear uncoupled springs may 
be adequate for real systems and it may be used to model 
the rolling friction phenomenon. 
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